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Abstract 

The calculation of interatomic magnetic exchange interactions entering the Heisen- 
berg model from the standpoint of the density functional theory (DFT) is outlined 
for two Fe-based molecular magnets: a trinuclear complex with a Schiff base lig- 
and, which makes an antiferromagnetically coupled frustrated system, and a model 
bipyrimidine-connected planar network of Fe ions. First-principles electronic struc- 
ture calculations are performed using the real-space method Siesta and the full- 
potential linearized augmented plane wave FLAPW method FLEUR, correspond- 
ingly. We discuss the application of fixed spin moment technique for preparing the 
system in a given magnetic configuration, and the effect of intraatomic Coulomb 
correlation, approximated by the LDA+U technique, on the values of interaction 
parameters. 
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1 Introduction 



Since the first systematizations of their chemistry and basic properties by, e.g., 
Kahn in 1993 [1], molecular magnets are now developing into a promising class 
of magnetic nanomaterials with a great potential in applications in areas such 
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as magnetic storage, quantum computing, or magneto-optical devices. An up- 
to-date review on microscopic properties (spin density, exchange interactions, 
magnetic anisotropy) of molecular magnets accessible in experiment and from 
first-principles calculations, is going to appear in Ref. [2]. Its major part is sum- 
marized as an internet publication [3] . In the present work we address the issue of 
extracting the interatomic exchange interaction parameters entering the Heisen- 
berg model form the density functional theory (DFT). We discuss two chemical 
systems of different complexity. The actual calculations are done by two differ- 
ent methods. Siesta [4] using compact and strictly confined atom-centered basis 
functions (see Ref. [5] for details), and the FLEUR [6] code, a reahzation of 
the highly accurate full-potential linearized augmented plane wave (FLAPW) 
method. Both calculations used the generalized gradient approximation to the 
exchange- correlation f unctionals . 

We skip the subtleties in the formulation of exchange parameters within the DFT; 
a detailed discussion can be found in Ref. [3]. Basically, we refer to total ener- 
gies in different specially prepared magnetic configurations, calculated from first 
principles. Sec. 2 covers the extraction of meaningful results for a Fe-trinuclear 
system by applying the fixed spin moment formalism. Sec. 3 addresses a model 
Fe-binuclear system and, specifically, the importance of intraatomic correlation, 
treated beyond the conventional DFT formalism. 



2 Fixed Spin Moment treatment of a Fe-trinuclecir system 

One of the ways to arrive at the values of interaction parameters in a magnetic 
system whose properties are well decribed by the Heisenberg model is to compare 
total energies for different magnetic configurations. This is rather straightforward 
if, say, ferromagnetic (FM) and antiferromagnetic (AFM) states are well defined 
and at least metastable. A more complicated case is a frustrated system, an 
example of which is made of three Fe atoms incorporated into organic molecule. 
We consider specifically Fe3(OAc)3L3, a synthesized by Boskovic et al. [7] derivate 
of the Schiff base H2L = salicylidene-2-ethanolamine. The molecular unit used in 
the electronic structure calculation by the Siesta method is shown in Fig. 2 (a) of 
Ref. [7] . The magnetization measurements reported in Ref. [7] indicate an (almost 
AFM) ground state with the total spin 5*= 1/2, coming about from the interaction 
of three s=5/2 spins of Fe(lll) ions. Experimentally, the interaction parameters 
J of the Heisenberg model, introduced in Ref. [7] as H = —2 JijSiSj, were 

Fc pairs 

reported to be necessarily different for an acceptable fitting of data, even as the 
molecule nearly maintains a trifold symmetry axis: Jij = —15.1 K, —13.6 K, and 
-12.4 K. 

We proposed earlier [8] to use the Fixed Spin Moment (FSM) method [9] for 
selecting a magnetic configuration of a molecular magnet (a six-center "ferric 
wheel" ) , improving simultaneously the stability of the electronic structure calcu- 
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Fig. 1. Results of the FSM calculation for [Fe3(OAc)3L3]. See main text and explana- 
tions [19] for details. 

lation. For the present system this idea becomes even more important, because it 
is no more obvious how to simulate a (nearly) AFM trial magnetic configuration 
in the DFT run, nor how to enforce an inversion of a selected local magnetic 
moment. On the contrary, by running a sequence of FSM calculations from the 
maximal nominal value of the total spin 5*= 15/2 (and higher) down to zero, one 
can scan the variety of all possible magnetic solutions within the DFT. 

The results of our FSM calculation are shown in Fig. 1. The total energy (left 
panel) has several local minima, the lowest one being for the FSM value of 5 
/i-B- The right panel shows that this situation corresponds to two Fe magnetic 
moments aligned in parallel and the third one aligned in the opposite direction. 
A competitive energy minimum is found for FSM = 15 /Xb, i.e., the ferromagnetic 
configuration. The solution with FSM = 1 /Xb has higher energy than the latter 
one by 74 meV, per Fe atom. It should be noted that, whereas the ground state 
of the system is experimentally found to have S=l/2, it must have a mixed 
nature over many quantum states, and hence does not correspond to any single- 
determinant DFT result. Nevertheless, the tendency for antiparallel coupling of 
Fe magnetic moments is correctly provided by our calculation; moreover the order 
of magnitude of the interaction constants can be correctly recovered: consistently 
with the above formula for the Heisenberg Hamiltonian, —Jk, (ii^l^TT _^TTi)^ (8s^) 
amounts to 57 K. The error (overestimation by a factor of about 4) in the DFT 
calculation is nearly the same as we found earlier for the "ferric wheel" [8]. A 
possible reason for such systematic discrepancy can be the underestimation of 
the intraatomic Coulomb correlation, to be discussed in the next Section. 

Two issues might need further discussion. The HOMO-LUMO gaps shown in the 
middle panel of Fig. 1 are different in the FSM formalism for two spin directions, 
because the fixing of the total spin amounts to an imposition of an external mag- 
netic field, moving apart the chemical potentials in two spin channels. Normally 
this would require the consideration of a Zeeman term when dealing with total 



3 



energies. However, for a system with a common band gap in both spin channels 
a unique chemical potential can be found for both spin directions, therefore the 
FSM scheme would simply fix one or another of the metastable magnetic con- 
figurations at no additional energy cost. One can see that, indeed, a common 
band gap can be found in all situations relevant for our discussion, i.e. with FSM 
= 1,5, and 15 /Xb, thus justifying the extraction of the J value above. Another 
observation concerns the local magnetic moments at Fe sites: as the FSM value 
is changing, and even as the saturation of magnetization (at FSM = 15 /^b) is 
achieved, the local Fe magnetic moments merely flip, essentially maintaining their 
magnitudes of 4 //b, and not 5 //b- The nominal value of s—5/2, associated with 
a single Fe atom, comes about due to magnetic polarization of ligands. This sit- 
uation is identical to that described earlier for "ferric wheels" [3,8,10]. In spite 
of the very different structure of organic ligands, both systems have similarities 
in the nearest neighbourhood and in the charge state of Fe atoms (Oe twisted 
octahedral coordination in "ferric wheels" vs. slightly distorted octahedral O5N 
coordination in the present system). 



3 Fe-binuclear system: Effect of intraatomic correlation 

We turn now to the discussion of the effect of on-site intraatomic Coulomb corre- 
lation on interatomic exchange parameters. Our objective was to find a compact 
and yet realistic system, for which an electronic structure calculation could be 
done with the use of an ultimately accurate method within the DFT, that is, the 
FLAPW method, and to apply a phenomenological correction by the LDA+U 
formalism [11] on top of it. Specifically, we apply the realization of the method 
in the FLEUR code [6]. We derive a simple model system from the Fe-binuclear 
complex described by Real et al. [12] being subject to previous studies [13]. We 
neglect for the moment the spin crossover property, which makes this system par- 
ticularly interesting, and make our calculation for the high-spin configuration on 
both Fe centers. The chemical formula of the substance is [Fe(bt)(NCS)2]2-bpym, 
with bt = 2,2'-bithazoline and bpym = 2,2'-bipyrimidine, and the molecular unit 
as shown on the left of Fig. 2. We also performed electronic structure calculation 
for this substance in its molecular and crystaUized form by the Siesta method, 
allowing structure relaxation, the results of this study will be reported elsewhere. 
For the present simulation, as the molecule is yet too large and too low-symmetric 
for an affordable FLAPW calculation, we simplify the molecular units, substitut- 
ing the bt fragments by bpym - that is realistic, as this is encountered in related 
compounds - see, e.g., [13]. Moreover we cut and connect the (NCS) tails into 
-N=C=N- chains that makes a dense 3-dimensional lattice out of planar Fe- 
bipyrimidine groups (Fig. 2, right panel). The nearest neighbourhood of each Fe 
center remains a slightly distorted Ng octahedron, as was the case in the original 
molecular unit. An idealization is that all N-Fe-N bonds make now right angles, 
and the repeated bipyrimidine fragments are coplanar. 
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Fig. 2. Molecular unit of the Fe-binuclear system Fe[(bt)(NCS)2]2bpym described in 
Ref. [12] (left) and a simplified model thereof, used in the calculation (right). The Fe 
atom is in the middle of a slightly distorted Ng octahedron. 

The calculated densities of states (per unit cell and per one Fe site) in Fig. 3 
clearly indicate a hybridization of FeSd with the 2p states of neighbouring nitro- 
gen atoms. Since the minority-spin Fe 3d states are not completely empty, the 
local magnetic moment is less than 5 /^b- As for the previous trinuclear system, 
the magnetic moment induced on nitrogen neighbours substantially contributes 
to the net value, associated to each Fe site. This resolves, at least in part, a 
controversy between the spin value s=5/2, expected at the Fe site from formal 
valence considerations, and the local magnetic moment of our calculation. The 
cumulative magnetic moment, distributed over the Fe site and its neighbours, 
behaves like a "rigid" spin in the sense of the Heisenberg model. However, one 
can expect that a stronger localization of this distributed magnetic moment at 
the Fe site will be enforced by an inclusion of intraatomic Coulomb correlation, 
underestimated in conventional DFT calculations. In order to check this possible 
effect, we brought in an additional Coulomb correlation applying a semi-empirical 




Energy (eV) 
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Fig. 3. Total density of states per unit cell (outer solid line) and local DOS per Fe site 
(evaluated within the muffin-tin sphere of 2.3 Bohr; shown by filled black) densities of 
states calculated for the model Fc-binuclear system, shown in the right panel of Fig. 2, 
for parallel orientation of Fe magnetic moments. Left panel: LDA calculation, right 
panel: LDA-|-?7 calculation with U =4 eV. 
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Table 1 

Local magnetic moments at the Fe site (in jj,^, within muffin-tin sphere of 2.3 Bohr); 
total magnetic moments per Fe center, total energy difference between FM and AFM 
configurations, and corresponding extracted values of the Heisenberg exchange param- 
eter J, assuming nominal spin value s=5/2. 





^Fe 


M/Fe 






U (cV) 


FM AFM 


FM 


- ^'afm (meV) 


J(K) 





3.62 3.61 


4.10 


102.5 


-95 


4 


3.93 3.92 


4.94 


76.8 


-71 



"LDA+C/" scheme [11]. We choose the Coulomb interaction constant U — 4: eV, 
which seems reasonable for correlated Fe compounds [14,15]. One sees in the 
DOS of Fig. 3 (right panel) the known effect of the "+f/" correction to lower the 
occupied states in energy whereas shifting unoccupied ones upwards. 

The values of magnetic moments and exchange parameters J, determined by the 
formula E = — JsiS2, are listed in the Table. It is noteworthy that, whereas 
the total magnetic moment per Fe site is increased in the LDA+t/ calculation, 
the local Fe magnetic moment still does not exceed 4 /Xb. On can conclude that 
the spatial distribution of magnetic density is not confined to the Fe atom even 
as the intraatomic Coulomb correlation is enhanced; instead, the magnetization 
maintains its complex spatial distribution, spilling onto neighbouring N atoms. 

As the intraatomic Coulomb correlation shifts apart on the energy scale the cen- 
ters of gravity of occupied and unoccupied states, it reduces the interatomic 
exchange parameter, because the latter (in one of its formulations, see [16]) has 
corresponding energy differences in the denominator. The same trend prevails in 
our case even as we estimate the values of J by comparing the total energies 
of two magnetic configurations. The effect of Coulomb correlation, included on 
top of a straightforward DFT calculation, is therefore always towards lowering 
the absolute values of Heisenberg exchange parameters. The same conclusion has 
been drawn by Boukhvalov et al. in their calculations for "Mni2" [17] and "V15" 
[18] molecular magnets. Comparing with the results of the previous Section, one 
can conclude that the estimation of exchange parameters in molecular magnets, 
that is already qualitatively correct (in sign and order of magnitude) in a conven- 
tional DFT calculation, can be brought closer to the experimental estimates by 
reducing their numerical values due to additionally included intraatomic Coulomb 
correlation. 
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